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Abstract. We study subdiffusive ratchet transport in periodically and randomly 
flashing potentials. Central Brownian particle is elastically coupled to surrounding 
auxiliary Brownian quasi-particles which account for the influence of viscoelastic 
\ environment. Similar to standard dynamical modeling of Brownian motion, the 

external force influences only the motion of central particle not affecting directly 

S the environmental degrees of freedom (see video). Just a handful of auxiliary 
Brownian particles suffice to model subdiffusion over many temporal decades. Time- 

■ modulation of the potential violates the symmetry of thermal detailed balance and 

induces anomalous subdiffusive current which exhibits a remarkable quality at low 
temperatures, as well as a number of other surprising features such as saturation at 
low temperatures, and multiple inversions of the transport direction upon a change of 
the driving frequency in nonadiabatic regime. Our study generalizes classical Brownian 

■ motors towards operating in sticky viscoelastic environments like cytosol of biological 
t"**- . cells or dense polymer solutions. 
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1. Introduction 

Physics of noise-assisted driven transport presents currently a well-established area 
of research pfl|2]. Most papers are devoted to classical transport in neglecting non- 
Markovian and even inertial effects. The corresponding stochastic nonlinear dynamics 
can be described by a Langevin equation in time-dependent potentials and/or by 
associated Fokker-Planck equation for the noise-averaged dynamics of the probability 
density of an ensemble of moving particles. Key ingredients are nonlinear dynamics 
in periodic potentials unbiased on average, friction and thermal noise related by the 
fluctuation-dissipation theorem (FDT), and an external driving which violates the 
symmetry of thermal detailed balance ensured by the FDT at thermal equilibrium. The 
emerged dissipative out-of-equilibrium directed transport is necessarily accompanied by 
an entropy production and related heat dissipation. 

Such profoundly out-of-equilibrium transport should be distinguished from other 
possibilities such as transport in a running potential, where the particle remains bound 
to a potential well which moves in space. Basically, in this later case one does not need 
even a periodic potential. Any trapping confining potential U (x) with a deep minimum 
at x = xq will convey a transport if this minimum is moving with velocity v, i.e. U(x) 
is replaced by U(x — vt). Such a "peristaltic" transport is clearly not related in essence 
to any entropy production and for a periodic running potential it can be symbolized 
by the Archimedean screw pump. This is a purely mechanical system, where both 
the friction and the noise are not principal for the transport occurrence. Quantum- 
mechanical counterparts of the Archimedean pump are also well-known [3] . A peristaltic 
pump can be also realized with highly overdamped isothermic thermodynamic systems 
operating close to thermal equilibrium, where dissipation keeps the particle fluctuating 
near to the potential minimum, and the energy lost due to friction is perpetually restored 
due to stochastic impact of the environment - the physical content of FDT. Friction 
plays here in fact a constructive role. Moreover, neglect of the energy gain from hot 
environment, at odds with FDT, is one of the common mistakes in the literature, leading 
to erroneous belief that one must always minimize dissipation and cool the environment 
in order to achieve at the highest efficiencies of isothermal engines possible. This is not 
necessarily so. For example, the so-called dissipationless Hamiltonian ratchets [I] cannot 
do any useful work at all against a load, i.e. are pseudo-ratchets with zero efficiency. A 
dissipative system can remain locally close to thermal equilibrium, but the form of the 
potential is cyclically and adiabatically slow changed by a driving force so that excess, 
uncompensated heat exchange between the particle and the thermal reservoir can be 
minimized. In this quasi-equilibrium scenario, the work required to change the form 
of potential (to move the minimum) can be transformed with minimal heat losses into 
the work against the load which opposes this potential modulation. This is why such 
an isothermal machine considered as a free energy transducer can in principle operate 
with the efficiency, defined as the ratio of the work against a load to the free energy 
spent to drive a working cycle, close to one [5HH]. For example, the efficiencies of highly 



Flashing subdiffusive ratchets 



3 




o o c# o o o 

c) 

o o o o 



d) 




Figure 1. Snapshots of numerical simulations of the dynamics in Eq. (fT3| at different 
instants of time a) t = 0.1, b) t = 0.4, c) t = 1.6, d) t = 2.0, e) t = 2.5. One Brownian 
particle (big filled circle) is coupled to auxiliary particles (small empty circles) at 
temperature T = 0.1 and moves in a periodically flashing with frequency v = 1.0 
potential ([T5]) with amplitude Uq = 0.75. 

optimized biological ionic pumps as high as 0.75 are common [9] and the efficiency of 
ATP synthase can approach the theoretical maximum of one [T0HT2] . 

The primary focus of this paper is different, on the thermal noise assisted transport, 
where the thermal noise is assumed to play a profound and constructive role. A 
paradigmatic example is provided here by transport in flashing potentials |13] . such 
as one in figure [TJ In a standard Markovian overdamped setup, when the potential is 
off, an initially localized particle diffuses with the position variance growing linearly, 
(5x 2 (t)) oc t, which corresponds to normal diffusion. When the potential is on, 
the particle relaxes to a minimum of the potential and the probability distribution 
becomes at least bimodal with a larger peak corresponding to sliding down a less steep 
side of the potential (as a larger basin of attraction corresponds to a larger distance 
between the minimum and maximum of a spatially asymmetric but periodic potential). 
When flashing repeats, either periodically or stochastically, but sufficiently slow, the 
net transport is expected to emerge in the left, natural direction, with the averaged 
particle's position growing linearly in time. Such normal transport is characterized by 
mean velocity. This natural direction is opposite to one in the fluctuating tilt potential 
ratchets of the same potential form [HIEE]- For a fixed unbiased on average potential, 
the total heat exchange between the particle and its environment is zero, FDT holds, and 
directed transport is forbidden by the symmetry of thermal detailed balance. Out-of- 
equilibrium potential fluctuations violate this symmetry and induce directed transport. 
There emerges an overall uncompensated excess heat flow to the environment associated 
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with the corresponding entropy production. A part of the energy put in the repeating 
potential flashing is dissipated as excess heat and a part can be used to do useful work 
against a load. If the load is absent all the consumed energy is dissipated as excess heat 
(futile motor) because the mechanical energy of the motor particle remains on average 
not changed. This is a well-established by now physical picture [16J. This basic model 
explains e.g. operating of single-headed kinesin motors [9J, where the energy to drive 
stochastic cycles is provided in effect by the energy of ATP hydrolysis. 

Generalization of this approach to account for a non-Markovian friction with 
memory is not trivial. The memory effects appear e.g. due to viscoelasticity of 
the environment [TTH23] . Even if a corresponding Generalized Langevin Equation 
(GLE) [24H27] is well known for any linear model of friction with the memory, the 
corresponding non-Markovian Fokker-Planck equation (NMFPE) [28--32J remains simply 
unknown for general nonlinear force-fields. The cases of constant force, or a linear in 
coordinate force, where the corresponding NMPPEs are known [28J, are not especially 
useful in the present context. Especially interesting is the case of a power law decaying 
memory kernel [18] which corresponds to subdiffusion and is associated with the Cole- 
Cole dielectric response of viscoelastic media [33j|3l]. Such viscoelastic subdiffusion has 
recently been found relevant also for transport in polymer networks [351 [36], cytosol 
of biological cells [37] and akin crowded fluids [38]. The corresponding fluctuating 
tilt or rocking ratchets in viscoelastic media have recently been proposed and studied 
in Refs. [SUED] with a number of quite unexpected and surprising properties revealed. 
The interest to ratchet effect in glass-like environments is growing [H] . Flashing ratchet 
transport in dynamically disordered potentials was proposed and studied earlier in 
Ref. [42]. Transport in disordered potentials is known to be equivalent within the 
mean field effective medium approximation to a continuos time random walk (CTRW) 
approach with independent residence times in traps [4l3H4"5] . This places our research 
into a general context of anomalous transport processes in complex media [4~3"]-46j . 
However, our approach to anomalous transport in viscoelastic media is very different 
from one based on uncorrelated CTRW. This work is devoted to flashing viscoelastic 
ratchets which turn out to be no less surprising than rocking viscoelastic ratchets [39jl4"0] 
and we shall take the advantage of the approach to anomalous subdiffusive transport 
developed recently in Refs. p2|[23l[39]. 

2. Dynamical approach to viscoelastic transport and stochastic modeling 

We start with a well-established dynamical approach to the theory of Brownian 
motion [251127] . The environment is modeled by a set of iVo particles with masses 
rrii harmonically coupled with the spring constants K{ to the Brownian particle of mass 
m: 




(1) 



(2) 
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where qi are the coordinates of the environmental particles, x is the coordinate of the 
Brownian particle, and f(x,t) is an external force which affects only the motion of 
Brownian particle but not the environment. As usually, using the Green function of 
harmonic oscillator one can express the bath oscillators coordinates <&(£) in Eq. (j2J) via 
the initial values qi(0) and Pi(0) for arbitrary x(t): 

Qi(t) — 9i(0) cosfcjjt) H — - — sinfwjt) H — / sinL>j(t — 

'/>,~\ Jo 

= [qM - x(0)] cos(wit) + sin(wit) + x(t) 

rriiUJi 

- / cos[w i (t-t')]ir(t / )tft . (3) 
Jo 



Here, = y Ki/rrii are the bath oscillators frequencies and the integration by parts has 
been used to obtain the second equality. By substituting fl3]) into ([I]) one immediately 
obtains 

-t 



where 



and 



mx+ [ r](t - t')x{t')dt' = f(x,t) +£(t), (4) 
Jo 

rj(t) = ^Kicosiuiit), (5) 



(6) 



Eq. (j3j) is still a purely dynamical equation of motion which is exact, the dynamics of 
irrelevant degrees of freedom is but excluded. For example, it describes a time-reversible 
dynamics, if external field f(x,t) does not violate the time-reversal symmetr)||]. Then 
the trajectories governed by Eq. (Tj0) are also most obviously time-reversal symmetric 
(see [17] on this point from stochastic perspective). This is contrary to a widespread 
misperception. The statistical irreversibility comes about on the level of a bunch 
of trajectories averaged over different initial realizations of qi{0), Pi(0). The loss of 
information due to averaging or coarse-graining is a well known source of irreversibility 
leading to statistical description of dynamical systems (see e.g. [19]). As a matter of 
fact, the dynamics described by Eq. is just a projection of a highly dimensional 
Hamiltonian dynamics onto two dimensional (x,p) phase subspace. 

Furthermore, one proceeds as in a typical molecular dynamics setup. The initial 
positions qi(0) and momenta Pi(0) of the environmental oscillators in Eq. ([6]) are sampled 
(this is the only non-dynamical element in the theory) from a canonical distribution at 

| Time-reversal symmetry can be dynamically broken by an external time-dependent field. For 
example, a harmonic mixing driving, f(t) = Ax cos(fM) + A 2 cos(2fM + 0) does violate the time-reversal 
symmetry dynamically for (f> ^ 48 . 
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temperature T, 



Mte(0),Pi(0)}|i(0)) = iexp 



1 E(^ + *.fe(0)-*(0)] 2 ) 

1=1 ^ 1 ' 



(7) 



2A; B T 

conditioned on the initial position of the Brownian particle x(0), where Z is the 
statistical sum of bath oscillators. Here one assumes that initial velocities of the bath 
oscillators are centered at zero, i.e. the medium is not moving as whole. Likewise, 
the average (<7i(0)) = x(0), i.e. the medium is initially equilibrated adjusting to the 
Brownian particle localized at x(0). The corresponding random force £(t) becomes a 
stochastic process which is obviously Gaussian and can be completely characterized 
by its first two statistical moments. First moment is obviously zero, (£(£)) = 
0. Furthermore, with Gaussian averages (pj(O)pj(O)) = SijiriiksT, (xj(O)xj(O)) = 
S^ksT/ni, (xi(0)pj(0)) = it is easy to show that 

mm) = k B Tr,(\t-t'\) (8) 

for any set of bath oscillators. This is the celebrated fluctuation-dissipation relation, 
or the second fluctuation dissipation theorem by Kubo [24J. The bath oscillators are 
conveniently characterized by the spectral density [2 



_ 2 



It allows to express r)(t) as r)(t) = {2 / tx) ^ du J {u) cosiuit) / uj and the noise 
spectral density as S(u) = 2k B TJ(u)/uj via the Wiener-Khinchin theorem, S(u>) = 
r^ o (^(t)^(0))e^*dt. The choice J(u) = r) a \ sin(vra/2)|u; a with < a < 2, yields the 
fractional Gaussian noise (fGn) £(£) introduced by Mandelbrot and van Ness [50]. For 
< a < 1 (sub-Ohmic thermal bath [27]). 

v(t) = V a t- a /T(l-a), (10) 

where T(z) is a standard gamma function. This choice yields subdiffusion 
asymptotically, {Sx 2 (t)) ~ t a , for an ensemble of particles [27]. The corresponding GLE 
PJ is termed also the fractional GLE, or FLE [231E21EI1E2] upon the use of the notion 
of fractional Caputo derivative to shorthand the frictional term, and the coefficient rj a 
is termed then the fractional friction coefficient. Such a dynamical modeling requires a 
very large iVo — > oo number of the bath oscillators with a quasi-dense spectrum. 



2.1. Stochastic modeling 

Alternatively, one can use just a handful of iV auxiliary Brownian particles clouding 
around the central particle [231ES] while modeling the rest as friction and noise acting 
on these representative ones. Such Brownian quasi-particles can serve to model sticky 
viscoelastic media. A particle clouded by other particles reminds conceptually polaron 
picture in condensed matter physics. Then, one replaces Eqs. (ITJ2"|) with 

N 

mi = f(x, t) - ^ kj(x - xi) (11) 
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rriiXi = ki{x- Xi) - 7]iXi + y/2rjik B TCi(t) , (12) 

where x^ are the coordinates of auxiliary particles, fcj are the corresponding coupling 
constants, ^/2rj i k B TC < i{t) are thermal Gaussian forces, with zero range of correlation, 
(Ci(t)(j{t')) = SijS(t — t'), and rji are the corresponding friction coefficients. Moreover, 
the overdamped limit for these auxiliary stochastic medium oscillators, rrii — > 0, yields 

x = v , 

N 

mv = f(x,t)+^2ui(t) (13) 

i=l 

itt = - kiV - ViUi + \/2v i k i k B T( i (t) , 

where v\ = ki/rji are the relaxation rates of the introduced viscoelastic forces u% = 
—ki(xi — x). The last equation for itj is similar to the Maxwell's relaxation equation 
for viscoelastic force in macroscopic theory of viscoelasticity [UJ which is augmented 
by the corresponding Langevin force in accordance with FDR. Such a description was 
introduced in Refs. [22l[23l|39l|53] to model anomalous Brownian motion in complex 
viscoelastic media within a generalized Maxwell model. For a particular case of the 
only one auxiliary particle (Maxwell model of viscoelasticity), the earlier description 
in Refs. [MJ|55] leading to GLE with thermal Ornstein-Uhlenbeck noise is readily 
reproduced. Excluding the dynamics of auxiliary variables Ui(t) and assuming that 
initially forces Wj(0) are thermally Gaussian-distributed with zero mean and dispersion 
(uf(0)) = kik B T yields again the GLE (Tj0) with the corresponding memory kernel 
presented by a sum of exponentials, 

N 

7?(t) = 5>e-^. (14) 

i=l 

The corresponding noise £(t) presents the sum of Ornstein-Uhlenbeck noises. 

The power-law memory kernel (fTOl) can be reliably approximated by such a sum [22] 
over a large time-interval [ti,th] by choosing z/j = v^jb 1 and ki = C a (b)r] a i / °' /T(l — a). 
Here, b > 1 is a dilation (scaling) parameter and C a (b) is a fitting dimensionless constant. 
Physical meaning of uq is a high-frequency (short-time) cutoff of the stochastic process 

ti = u^ 1 . The low-frequency (long-time) cutoff corresponds to th = b N ~ 1 t i . Similar 
scaling and approximation to a power law are well known in the theory of anomalous 
relaxation [13],[56]. By adjusting b and iV one can approximate power law over about 
r = N log 10 6—2 time decades. Weak dependence of r on b and iV ensures a very powerful 
and accurate numerical approach to FLE dynamics. Just a handful of auxiliary particles 
can suffice for all practical purposes. Of course, such a Markovian embedding of FLE 
is not unique [57] . However, it appeals by a clear physical interpretation. Recently, the 
set ([TBI was formally generalized to numerically integrate also superdiffusive FLE [58J. 
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3. Flashing ratchet 

To study flashing ratchet transport we take f(x,t) = fo(x)r(t). Here, fo(x) = —U'{x) 
is a deterministic force acting on the particle (Brownian motor) when a potential U(x) 
is switched on. Furthermore, r(t) switches on/off taking just two values, zero and one, 
r(t) = {1,0}. For U(x) we take a typical ratchet potential 



with amplitude Uq and period L. For r(t) two models are considered: periodic vs. 
random switching. In the periodic case, r(t) = 1/2 [1 + sign (sin(27rz/t))], where sign(-) 
is standard signum function and v is linear frequency of oscillations. The first switching 
occurs at t\ = l/2u and the total number of switching is N sw = T tot jt\ = 2uT tot , where 
T tot is a total computation time which equals N s dt; N s is a total number of integration 
steps of Eq.fll]), dt is a time-step. Finally, one can write N sw = 2vN s dt. For random 
switching r(t), we consider the model of dichotomous Markovian process (DMP) with 
the probability p = 2vdt to switch within dt. Such DMP is symmetric with equal 
residence time distributions in two states, ip(r) = 2v exp(— 2vt). Here, 2v is transition 
rate from one state to another one (and also inverse of the mean residence time in a 
state). The mean flashing frequency (v) equals v. 

In all simulations we scale coordinate in units of L and time in units of t v = 
(m/r] a ) 1 ^ 2 ^ a \ It is assumed to be temperature- independent in accordance with the 
underlying Hamiltonian model j25j[27]. This is a standard assumption done also in 
other ratchet models [UI2]. Furthermore, energy is scaled in E = L 2 r]a^ 2 ~ a ^ /m a ^ 2 ~ a ' ) 
and temperature T = E/ks- In this work we choose a = 0.5 and b = 10, Ci/2(10) = 1.3, 
and use z/ = 100, iV = 12. Such a set of parameters gives an excellent approximation to 
the FLE dynamics over at least ten time decades, until t max = 10 8 . Similar to [221139] this 
was checked by comparison with the exact analytical solution for the position variance 
obtained within GLE and FLE [2ljl5T| in the force-free case. The numerical errors due to 
the memory kernel approximation are negligible as compare with the overall statistical 
error in our stochastic simulations (several percents), see for details in [22]. Like in |40j, 
numerical solutions of stochastic differential equations were performed on the graphical 
processor units (GPUs) with double precision using stochastic Heun method. The use 
of GPU computing allowed to parallelize and accelerate simulations by a factor of about 
100, as compare with conventional CPU computing. 

Schematically, viscoelastic dynamics in the flashing ratchet potential is shown in 
figured] (see also video abstract). The auxiliary particles are placed on the level U = 
and they are not influenced by potential. Initially, when the potential is switched on, 
the Brownian particle moves in the potential well (see snapshots a and b). When the 
potential is off (snapshots c and d) the Brownian particle "freely" subdiffuses, interacting 
only with auxiliary particles. After the potential is switched on again, Brownian particle 
moves into the next potential well (see snapshot e) with the net motion occurring in the 
negative direction. During such motion not all auxiliary particles follow immediately to 




(15) 
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the Brownian particle. Some of them are at large distances from the Brownian particle 
and are essentially less mobile than the Brownian particle and its nearest environment 
because they have essentially larger friction coefficient than other auxiliary particles. 
The weaker is the coupling constant fc, of an auxiliary particle to the Brownian one the 
larger is the corresponding frictional coefficient r/j. These very sluggish particles create 
a slowly fluctuating quasi-elastic force (a temporal biasing force [22],[23]) acting on the 
central particle. 

4. Results and discussion 

In all simulations we monitored two main statistical quantities, the mean position 
(x(t)) and the position variance (d~x 2 (t)) = (x 2 (t)) — (x(t)) 2 . For the ensemble average, 
iVens = 10 4 trajectories were used. The focus is on such physical quantities as subdiffusive 
current (subvelocity) v a , subdiffusion coefficient D a and generalized Peclet number [39] 
Pe Q := v a /D a . The latter one surves as a measure for the coherence and quality of 
anomalous stochastic transport, by analogy with normal one [61)]. The subvelocity and 
the subdiffusion coefficient are defined as usually [23 | [39 l l43 | l53] 



Here, the limit is understood in the following sense: t is large but still much smaller than 
the memory cutoff, which can always be made unreachable in numerical simulations and 
thus is irrelevant for the results presented. Practically, the values of v a and D a were 
calculated by fitting the numerical dependencies (x(t)) and (Sx 2 (t)) with a power-law 
function a^ v j^yt a , extracting the corresponding a v and within the last time window of 
simulations. Total time of simulations was varied in the interval [2 x 10 5 , 10 6 ] depending 
on system parameters to guarantee a good fit with a = 0.5 (convergency is slow); the 
time step was fixed at dt = 2 x 10~ 3 . 

Typical dependencies for the mean particle position and variance on time for some 
different values of the potential amplitude U and flashing frequency v are shown in figure 
E|a) and (b), respectively. Here black solid, red dash and blue dash-dot lines correspond 
to the periodic flashing, whereas green dash-dot-dot line relates to the random one. The 
direction of transport in figure [2(a) is negative and the transport has clearly subdiffusive 
character. We indicate in figure [2](a) also the corresponding power- law asymptotics. One 
can notice that stochastic flashing tends to delay the corresponding transport process 
in comparison with the periodic one of the same frequency (cf. red dash vs. green dash- 
dot-dot lines), and the transport becomes faster with the increased flashing frequency 
(this is but not a universal feature, see below). The diffusion is initially always ballistic 
(universal regime), (5x 2 (t)) = v^t 2 , where vt = ^/ksT '/m is thermal velocity. This is 
because the Brownian particle's velocity is initially thermally distributed and the action 
of the medium requires some time to settle in. After some transient, subdiffusion follows 




(16) 



(17) 




Figure 2. The mean particle position (a) and mean variance (b) at different values 
of the potential height Uq and flashing frequency v for periodic and random (green 
dash-dot-dot line) driving at T = 0.25. The black line in (b) depicts exact analytical 
FLE result for the potential-free case. 



asymptotically to one and the same universal dependence, (Sx 2 (t)) = 2D a t a /T(l + a) 
with D a = kBT/r] a , as in the absence of potential, independently of the potential 
height and the presence of driving, cf. in figure E^b). This universality (or a weak 
sensitivity in the case of a strong fast driving [39]) is a benchmark of viscoelastic 
subdiffusion j22j|23l[39l|53] . We elaborate on this fact in more detail below. 

4-1. Ergodicity 

An important issue in anomalous transport is ergodicity (EH [62], i.e. whether time- 
average over a single particle trajectory delivers non-random result, same for all identical 
particles, or a principal randomness and scatter in the single-particle averages emerge. 
If this is the case, then the behavior and fate of each individual though identical particle 
is different, even in the limit of infinitely long trajectories, and only ensemble-averaging 
smears out this principal randomness and delivers non-random result. In the case of 
ergodic transport, a single-trajectory average should coincide with the result of the 
ensemble averaging. For example, CTRW subdiffusion with divergent mean residence 
times in traps is patently nonergodic [63] and single-trajectory averages obey some 
universal fluctuation laws [63H66] leading to a universal scaling relating such transport 
and diffusion in periodic potentials [68]. On the contrary, viscoelastic subdiffusion, 
both free [67], and in time-independent potentials [22] was shown to be mostly ergodic, 
though some transient nonergodic features can be present. Whether it remains ergodic 
also in time- fluctuating fields presents a nontrivial problem. First, we checked the 
asymptotic ergodicity of ratchet viscoelastic transport which is expected. To show that 
this indeed is the case we plotted the fluctuating subvelocity v a (t), obtained for a single 
trajectory as v a {t) = T(l + a)x(t)/t a in figure [3](a) vs. the ensemble-average obtained 
with 10 4 particles. One can see that v a (t) fluctuates around the ensemble average 
with the amplitude of fluctuations diminishing with t. In the limit t — > oo the both 
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a) i h) t 

Figure 3. Subvelocity v a (t), obtained for the single trajectory a) for periodic and 
random flashing ratchets at Uq — 0.5, T — 0.25 and linear frequency of flashing (its 
mean value) v — 0.5. Mean deviation (5x 2 (t)} b) for ensemble (E) and single trajectory 
(S) averaging at Uo = 0.5 and T = 0.25 for periodic and random flashing ratchets. 
The time averaging for a single trajectory is done using 16-exponential approximation 
to memory kernel and T = 2 x 10 7 . 



averages coincide and the transport is ergodic. For any finite t, statistical fluctuations 
are present, of course, also in ergodic case. Next, a single trajectory average of the 
squared displacement can be obtained as [22 | l6l ^ l65 | l6T]. 

(Sx 2 (t)) T = £ W + - x(t')] 2 dt>. (18) 

In the ergodic case, this average should coincide in the limit T — > oo but for any t 
with the result of the ensemble averaging. If the agreement holds only for large t, 
then diffusion is asymptotically ergodic. For viscoelastic subdiffusion in fixed periodic 
potentials the agreement holds both for small t (ballistic regime) and for large t [22J. 
However, some deviations from ergodicity occur for intermediate t, on the time scale of 
diffusion over several potential periods. This is because even if the mean escape time 
exists the escape kinetics remains anomalous for activation barriers of an intermediate 
height [22]- 

Diffusional spreading derived from single trajectories is compared with the 
corresponding ensemble averages in figure [S^b). In the ballistic regime, using x(t + 
t') — x(t') ~ v(t')t for small t in Eq. ( |T8|) . one obtains 

(5x 2 (t)) T ^t 2 ^— t £ t v\t)dt= (Sv 2 (t)) T t 2 . (19) 

Clearly, if in the limit T — > oo the time average of v 2 (t) coincides with the ensemble 
average yielding v^, then ergodicity holds also in the ballistic regime. In the absence 
of driving this is indeed the case [22]. However, a periodic driving exciting coherent 
oscillations emerging due to a combined action of the external trapping potential and 
the viscoelastic cage effect can cause large additional temporal fluctuations in v(t) which 
are not self-averaged in {v 2 (t))-r. Accordingly, such a periodically driven subdiffusion 
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Figure 4. Anomalous current (subvelocity v a ) as a function of linear frequency of 
periodic flashing at fixed temperature T = 0.25 for different values Uo. 

ceases to be ergodic in the ballistic regime. Temporally but on a large intermediate time 
scale ergodicity is generally broken. Nevertheless, asymptotically it remains ergodic as 
figure [3^b) implies. Convergency to such asymptotically ergodic regime is very slow. 
However, it improves with the increased frequency of flashing, see the case v = 5 in 
figure [3](b) . For random driving, ergodicity in the discussed sense holds also in ballistic 
regime. 

4-2. Frequency dependence 

We proceed further with the dependence of anomalous current v a on frequency v at 
a fixed value of temperature for different potential amplitudes Uo- In the case of 
fluctuating tilt subdiffusive ratchets (39] this dependence was especially intriguing. That 
anomalous transport vanished in the limit of vanishing modulation frequency (adiabatic 
driving limit) and exhibited a maximum at intermediate frequencies. The origin of 
this maximum was related to stochastic resonance effect in Ref. [40J. This is in a 
sharp contrast with normal diffusion fluctuating tilt ratchets where transport optimizes 
namely in the adiabatic limit. The results for periodic flashing are shown in figure 
H] and for stochastic driving in figure [5j For small driving frequencies subdiffusive 
transport occurs always in the natural negative direction, which is opposite to the 
natural positive direction of fluctuating tilt ratchet. Our anomalous Brownian motors 
share these features with their normal diffusion counterpart [HE]- There are but novel 
features which are due to a profound role of the inertial effects in our model. We recall 
that in the used scaling the unit of time corresponds to a scaling time constant t v of 
the decay of velocity autocorrelation function. Correspondingly, the ballistic regime 
occurs until t ~ 1, when the potential is off. For periodic driving and a small potential 
height Uq = 0.25 = T (see red line with squares in figure @J subvelocity takes negative 
values for all flashing frequencies considered. The dependence of \v a \ on v has a unique 
maximum at z/ c q ~ 0.5, where subdiffusive transport is optimized, with negative v a . For 



Flashing subdiffusive ratchets 13 





-0.02 



-0.04 



-0.06 

0.5 1.0 1.5 2.0 2.5 

<V> 

Figure 5. Anomalous current (subvelocity v a ) as a function of mean linear frequency 
of random flashing at fixed temperature T — 0.25 for different values Uq. 

the fast-flashing regime, when v > 1 subvelocity asymptotically decays to zero with v. 
Increase in the ratchet potential height 2Uq accompanied by an increased role of the 
inertial effects leads to a profound change in the character of the frequency-dependence 
Vaiy) (see black curve with circles for Uq = 0.5 in figure H]). Let us study this dependence 
in more detail. The first difference from the former case is a multiple (double) current 
inversion, realized around z/ c i ~ 0.335 and u C 2 ~ 0.402. Enlarged picture of this inversion 
effect is shown in the insert a) placed on the right panel in figure HI Here, a well defined 
peak occurs at z/ max i = 0.375 with positive values for v c \ <v < u c2 . Furthermore, on the 
left side of v < z/ cl there is a maximum of \v a \ for negative v a . Increase of v beyond z/ c2 
leads to a global maximum of \v a \ realized at ^ max 2 — 0.75. With a further increase in 
flashing frequency a complicated oscillatory pattern with several more relative minima 
and maxima in |t> a | emerges. With a further increase of the potential amplitude these 
oscillatory resonance-like features become more pronounced - compare two cases with 
Uq = 0.50 and Uq = 0.75 in figure HI This is clearly related to an increased role of the 
inertial effects. First, the amplitude of oscillations in v a {v) increases. Second, ever more 
extrema appear, see insets b) and c) in figure HI 

The situation is quite different in the case of random flashing. Dependencies of 
anomalous current versus mean linear frequency (u) at fixed temperature T = 0.25 
for different values of C/ are shown in figure |5j Here, contrary to the case of periodic 
flashing one has more simple dependencies v a {{v)) with one broad minimum for each 
considered potential amplitude Uq. Increase in Uq leads to increase in absolute values 
for subvelocity, broadening the minimum and shifting minimum position towards larger 
values of mean linear frequency Thus, in the random flashing case the system is not 
manifestly sensitive to inertial effects due to stochastic nature of flashing. The maximal 
value of subvelocity becomes also much smaller, compare figure H] and figure EJ due to 
the absence of resonances. 

A benchmark feature of viscoelastic subdiffusion is that it asymptotically does 
not depend on the presence of static periodic potential [22i[23j|53]. It was also only 
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Figure 6. Sub-diffusion coefficient D a as a function of linear frequency v (or 
mean linear frequency (v)) for periodic and random flashing, respectively, at fixed 
temperature T — 0.25 and different values Uq. 



weakly dependent on driving in the case of rocking subdiffusive ratchets [3H]- The 
diffusional spread in flashing ratchets is also practically not dependent on the potential 
amplitude Uq and flashing frequency as figure El illustrates. One can see, that for both 
periodic and random flashing the anomalous diffusion coefficient does not display any 
profound dependence on the flashing frequency and the potential amplitude. It takes 
values around 0.25, which is D a = T in the scaled units, with a deviation which is less 
than 3% and lies within the statistical error margin of our simulations. Therefore, the 
corresponding generalized Peclet number Pec, := v Q /D a , which measures the coherence 
and quality of transport, resembles the behavior of the absolute value of subcurrent in 
figures H] and El which should be multiplied by factor 4 in order to obtain Pe a . 

4-3. Temperature dependence 

The temperature dependence of subvelocity v a , and generalized Peclet number Pe a is 
mostly intriguing, whereas the dependence of subdiffusion coefficient D a on temperature 
is expected to be trivial, D a (T) = T. We plotted numerical results for all these 
three quantities in figure [7] (periodic flashing) and figure M (stochastic flashing) for 
three different test values of the flashing frequency, v = 0.1 (low- frequency), v = 
0.375 (intermediate frequency), and v = 0.75 (high-frequency). The temperature 
dependencies of subvelocity share two common features, see figures [7(a) andET.a). First, 
the subvelocity is finite in the limit T — > (this is always so for random flashing, but only 
for certain frequency windows in the case of periodic flashing). Second, it diminishes to 
zero in the limit of high temperatures. The latter is expected, but the former is not, being 
highly surprising. Indeed, in a combination with the subdiffusion coefficient following to 
the expected linear dependence and vanishing in the limit T — > [see figures [7Jb) and 
|8(b)], the non-zero value of v a in this limit means that the basic mechanism of operating 
our subdiffusive viscoelastic Brownian motors is very different from one in the case of 





Figure 7. Anomalous current (subvelocity) v a (a), subdiffusion coefficient D a (b) 
and generalized Peclet number Pe Q (c) versus temperature for the periodic flashing 
ratchets at different values for linear frequency of flashing v and potential height Uq. 
The legend for (a), (b) and (c) figures is the same. 



normal Brownian motion. This is at odds with our expectations (see Introduction). 
Indeed, in the absence of viscoelastic effects the flashing Brownian motors of the studied 
kind are known to be essentially based on thermal diffusion, i.e. such transport is assisted 
by thermal fluctuations [9]. It is therefore expected to vanish in the limit T — > 0, 
where the thermal fluctuations vanish (in neglecting quantum effects). Our results 
say, however, that viscoelastic flashing ratchets are not based primarily on the thermal 
diffusion and thermal fluctuations. These are long-range elastic correlations which are 
at work here in combination with the potential fluctuations. Namely, the external field 
acting on the Brownian particle changes abruptly (on-off process) the energy stored 
in the elastic springs. Not only the Brownian particle moves under the influence of 
fo(x), but also the auxiliary particles under the influence of Brownian particle. Some 
of them rapidly adjust their positions, but there are also slow particles which cannot 
follow immediately. When the potential is off these extra sluggish particles pull the 
Brownian particle. This is why its motion is generally not frozen even at T = within 
our purely classical setup. Because subdiffusion diminishes with lowering temperature, 
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Figure 8. Anomalous current (subvelocity) v a (a), subdiffusion coefficient D a (b) and 
generalized Peclet number Pe Q (c) versus temperature for the random flashing ratchets 
at different values for mean linear frequency of flashing (v) and potential height Uq. 
The legend for (a), (b) and (c) figures is the same. 



the generalized Peclet number grows accordingly, Pe Q (T) ~ 1/T 5 with 5 in the range 
[0.83, 1.05], see figure [7(c), and figure [8(c). Clearly, Pe Q — > oo in the limit T — > 0, where 
the transport becomes perfect. 

There are also profound differences in the temperature-dependence of anomalous 
transport in the cases of periodic and random flashing. So, in the case of periodic 
flashing the transport can occur in the counterintuitive positive direction, see the curve 
for Uq = 0.5, v = 0.375 in figure [7(a). In this case, v a remains finite also at T = 0, 
and an inversion of the current direction with increasing temperature is possible. This 
is not so for other three curves in figure [7(a), where transport vanishes jump-like at 
T = 0. Therefore, it seems obvious that for a periodic flashing there are frequency 
windows for permitted transport at T = (a detailed investigation of this issue is left 
for a separate study). On the contrary, for random driving in figure a, the transport 
remained always finite at T = 0. Moreover, here transport always occurs in the intuitive 
negative direction. Both in the random case and in the periodic case the transport can 
be optimized at some temperature T max 7= 0. It can be but also that T max = 0. 
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5. Conclusions 

With this work we put forward anomalously slow Brownian motors based on flushing 
subdiffusion in viscoelastic media. Both anomalous transport and subdiffusion were 
shown to be asymptotically ergodic. However, such driven subdiffusion can be 
transiently nonergodic on appreciably long time scales. The transport exhibits a 
remarkably good quality at low temperatures. It remains finite even at zero temperature 
in many cases (but not always in the case of periodic flashing) with zero dispersion. 
In this limit, the transport is clearly induced by long-range elastic correlations in a 
combination with potential flashing and not by the thermal noise of environment. This is 
surprising and at odds with a popular explanation of the origin of noise-induced flashing 
transport in the case memory-free stochastic dynamics [HI2]. Moreover, in the case of 
periodic driving anomalous current exhibits multiple resonance-like features, can flow in 
the counter-intuitive direction, and invert its direction both with the change of flashing 
frequency and with change of temperature. The inertial cage effects are essential for 
many observed features. In particular, they are crucial for the resonance-like features 
manifested for high potential barriers in the case of periodic flashing. We hope that 
our research will stimulate further cross-fertilization of the field of fluctuation-induced 
transport and the field of anomalously slow transport and subdiffusion which flourish 
at present with little interaction. There emerges an increasing experimental support for 
the occurrence of subdiffusion in cytoplasm of biological cells and we believe that there 
might be also place for operating sluggish molecular motors of the kind described. 
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